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Abstract. The pscudorapidity density and multiplicity distribution of charged particles produced in proton- 
proton collisions at the LHC, at a centre-of-mass energy ^/s = 7 TeV, were measured in the central 
pseudorapidity region |?7| < 1. Comparisons are made with previous measurements at ^/s = 0.9 TeV 
and 2.36 TeV. At = 7 TeV, for events with at least one charged particle in \ri\ < 1, we obtain 
dAfch/d?? = 6.01±0.01{stat.)toA2isyst.). This corresponds to an increase oi 57.6 %±0. 4 %{stat.)tf,l %{syst.) 
relative to collisions at 0.9 TeV, significantly higher than calculations from commonly used models. The 
multiplicity distribution at 7 TeV is described fairly well by the negative binomial distribution. 
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Introduction 

We present the pseudorapidity density and the muhiphc- 
ity distribution for primary charged particle^ from a sam- 
ple of 3 X 10^ proton-proton events at a centre-of-mass en- 
ergy y/s = 1 TeV collected with the ALICE detector [T] at 
the LHC [2] , and compare them with our previous results 
at = 0.9 TeV and = 2.36 TeV OS]. The present 
study is for the central pseudorapidity region jr/j < 1. 

In the previous measurements, the main contribution 
to systematic uncertainties came from the limited knowl- 
edge of cross sections and kinematics of diffractive pro- 
cesses. At 7 TeV, there is no experimental information 
available about these processes; therefore, we do not at- 
tempt to normalize our results to the classes of events 
used in our previous publications (inelastic events and 
non-single-diffractive events). Instead, we chose an event 
class requiring at least one charged particle in the pseu- 
dorapidity interval |?7| < 1 (INEL>0|^|<i), minimizing the 
model dependence of the corrections. We re-analyzed the 
data already published at 0.9 TeV and 2.36 TeV in or- 
der to normalize the results to this event class. These 
measurements have been compared to calculations with 
several commonly used models [SHlOj which will allow a 

^ Primary particles are defined as prompt particles produced 
in the collision and all decay products, except products from 
weak decays of strange particles. 



better tuning to accurately simulate minimum-bias and 
underlying-event effects. Currently, the expectations for 
7 TeV differ significantly from one another, both for the 
average multiplicity and for the multiplicity distribution 
(see e.g. [Tl]l. 

ALICE detector and data collection 

The ALICE detector is described in [T]. This analysis 
uses data from the Silicon Pixel Detector (SPD) and the 
VZERO counters, as described in [3l|4]. The SPD detector 
comprises two cylindrical layers (radii 3.9 cm and 7.6 cm) 
surrounding the central beam pipe, and covers the pseu- 
dorapidity ranges |?7| < 2 and \ri\ < 1.4, for the inner and 
outer layers, respectively. The two VZERO scintillator ho- 
doscopes are placed on either side of the interaction region 
at z = 3.3 m and z — —0.9 m, covering the pseudorapidity 
regions 2.8 < rj < 5.1 and —3.7 < rj < —1.7, respectively. 

Data were collected at a magnetic field of 0.5 T. The 
typical bunch intensity for collisions at 7 TeV was 
1.5 X 10^*^ protons resulting in a luminosity around 
10^'' cm~^s~^. There was only one bunch per beam collid- 
ing at the ALICE interaction point. The probability that 
a recorded event contains more than one collision was es- 
timated to be around 2 x 10"'^. A consistent value was 
measured by counting the events where more than one 
distinct vertex could be reconstructed. We checked that 
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pileup events did not introduce a significant bias using a 
simulation. 

The data at 0.9 TeV and 7 TeV were collected with 
a trigger requiring a hit in the SPD or in either one of 
the VZERO counters; i.e. essentially at least one charged 
particle anywhere in the 8 units of pseudorapidity. At 
2.36 TeV, the VZERO detector was turned off; the trig- 
ger required at least one hit in the SPD {\ri\ < 2). The 
events were in coincidence with signals from two beam 
pick-up counters, one on each side of the interaction re- 
gion, indicating the passage of proton bunches. Control 
triggers taken (with the exception of the 2.36 TeV data) 
for various combinations of beam and empty-beam buck- 
ets were used to measure beam-induced and accidental 
backgrounds. Most backgrounds were removed as described 
in (31. The remaining background in the sample is of the 
order of lO"'^ to 10"^ and can be neglected. 



Event selection and analysis 

The position of the interaction vertex was reconstructed 
by correlating hits in the two silicon-pixel layers. The ver- 
tex resolution achieved depends on the track multiplicity, 
and is typically 0.1-0.3 mm in the longitudinal (z) and 
0.2-0.5 mm in the transverse direction. 

The analysis is based on using hits in the two SPD 
layers to form short track segments, called tracklets. A 
tracklet is defined by a hit combination, one hit in the 
inner and one in the outer SPD layer, pointing to the 
reconstructed vertex. The tracklet algorithm is described 
in [3111]. 

Events used in the analysis were required to have a 
reconstructed vertex and at least one SPD tracklet with 
1 77 1 < 1. We restrict the z-vertex range to \z\ < 5.5 cm to 
ensure that the 77-interval is entirely within the SPD ac- 
ceptance. After this selection, 47000, 35 000, and 240 000 
events remain for analysis, at 0.9 TeV, 2.36 TeV, and 
7 TeV, respectively. The selection efficiency was studied 
using two different Monte Carlo event generators, PYTHIA 
6.4.21 |5i6j tune Perugia-0 [9, and PHOJET with de- 
tector simulation and reconstruction. 

The number of primary charged particles is estimated 
by counting the number of SPD tracklets, corrected for: 

— geometrical acceptance and detector and reconstruc- 
tion efficiencies; 

— contamination from weak-decay products of strange 
particles, gamma conversions, and secondary interac- 
tions; 

— undetected particles below the 50 MeV/c transverse- 
momentum cut-off, imposed by absorption in the ma- 
terial; 

— combinatorial background in tracklet reconstruction. 

The total number of collisions corresponding to our data 
is obtained from the number of events selected for the 
analysis, applying corrections for trigger and selection ef- 
ficiencies. This leads to overall corrections of 7.8 %, 7.2 %, 
and 5.7% at 0.9 TeV, 2.36 TeV, and 7 TeV, respectively. 



The multiplicity distributions were measured for I77I < 
1 at each energy. For the 0.9 TeV and 2.36 TeV data we did 
not repeat the multiplicity-distribution analysis, we use 
the results from [3] while removing the zero-multiplicity 
bin. At 7 TeV, we used the same method as described 
in [3J[T3] to correct the raw measured distributions for ef- 
ficiency, acceptance, and other detector effects, which is 
based on unfolding using a detector response matrix from 
Monte Carlo simulations. The unfolding procedure applies 

minimization with regularization |12j . Consistent re- 
sults were found when changing the regularization term 
and the convergence criteria within reasonable limits, and 
when using a different unfolding method based on Bayes's 
theorem [TilfTS] . 

Systematic uncertainties 

Only events with at least one tracklet in \ri\ < 1 have 
been selected for analysis in order to reduce sensitivity 
to model-dependent corrections. However, a fraction of 
diffractive reactions also falls into this event category and 
influences the correction factors at low multiplicities. In 
order to evaluate this effect, we varied the fractions of 
single-diffractive and double-diffractive events produced 
by the event generators by ±50 % of their nominal values 
at 7 TeV, and for the other energies we used the varia- 
tions described in [3]. The resulting contributions to the 
systematic uncertainties are estimated to be 0.5%, 0.3%, 
and 1 % for the data at 0.9 TeV, 2.36 TeV, and 7 TeV, 
respectively. For the same reason, the event selection effi- 
ciency is sensitive to the differences between models used 
to calculate this correction. Therefore, we used the two 
models which have the largest difference in their multi- 
plicity distributions at very low multiplicities (see below) : 
PYTHIA tune Perugia-0 and PHOJET. The first one was 
used to calculate the central values for all our results, and 
the second for asymmetric systematic uncertainties. The 
values obtained for this contribution are -fO.8%, -1-1.5%, 
and -f2.8% for the three energies considered. 

Other sources of systematic uncertainties, e.g. the par- 
ticle composition, the px spectrum and the detector effi- 
ciency, are described in [3], and their contributions were 
estimated in the same way. As a consequence of the smaller 
uncertainties on the event selection corrections the total 
systematic uncertainties are significantly smaller than in 
our previous analyses, which use as normalization inelastic 
and non-single-diffractive collisions. Many of the system- 
atic uncertainties cancel when the ratios between the dif- 
ferent energies are calculated, in particular the dominating 
ones, such as the detector efficiency and the event gener- 
ator dependence. The systematic uncertainty related to 
diffractive cross sections was assumed to be uncorrelated 
between energies. 

Results 

The pseudorapidity density of primary charged particles 
in the central pseudorapidity region jr/j < 1 are presented 
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Table 1. Charged-particle pseudorapidity densities at central pseudorapidity (I77I < 1), for inelastic collisions having at least 
one charged particle in the same region (INEL>0|^|<i), at three centre-of-mass energies. For ALICE, the first uncertainty is 
statistical and the second is systematic. The relative increases between the 0.9 TeV and 2.36 TeV data, and between the 0.9 TeV 
and 7 TeV data, are given in percentages. The experimental measurements are compared to the predictions from models. For 
PYTHIA the tune versions are given in parentheses. The correspondence is as follows: D6T tune (109), ATLAS-CSC tune (306), 
and Perugia-0 tune (320). 



Energy 


ALICE 


PYTHIA [5l[6] PHOJET [TO] 


(TeV) 




(109) 7 (306) ,8 (320) [9] 




Charged-particle pseudorapidity density 


0.9 
2.36 
7 


3.81 ± O.Ollg g^ 
4.70 ± O.Oll°;Ji 
6.01 ±0.01^^;^^ 


3.05 3.92 3.18 
3.58 4.61 3.72 
4.37 5.78 4.55 


3.73 
4.31 
4.98 


Relative increase (%) 


0.9-2.36 
0.9-7 


23.3±0.4tJ:7 
57.6±0.4l?;^ 


17.3 17.6 17.3 
43.0 47.6 43.3 


15.4 
33.4 



PHOJET 

PYTHIA 

Pemgia-O 

PYTHIA 

ATLAS-CSC 

PYTHIA 

D6T 



ALICE 


INEL>0 |ri| < 1 




0.9 TeV ^2.36 TeV 


■ 


0.9 TeV ^7.0 TeV 



I 



20 



40 



60 

Increase (%) 



Fig. 1. Relative increase of the charged-particle pseudorapid- 
ity density, for inelastic collisions having at least one charged 
particle in \ri\ < 1, between ^ = 0.9 TeV and 2.36 TeV (open 
squares) and between = 0.9 TeV and 7 TeV (full squares), 
for various models. Corresponding ALICE measurements are 
shown with vertical dashed and solid lines; the width of shaded 
bands correspond to the statistical and systematic uncertain- 
ties added in quadrature. 
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Fig. 2. Charged-particle pseudorapidity density in the cen- 
tral pseudorapidity region |?7| < 0.5 for inelastic and non- 
single-diffractive collisions [4l ll6tf25] . and in \ri\ < 1 for inelas- 
tic collisions with at least one charged particle in that region 
(INEL>0|^l<i), as a function of the centre-of-mass energy. The 
lines indicate the fit using a power-law dependence on energy. 
Note that data points at the same energy have been slightly 
shifted horizontally for visibility. 



in Table [T] and compared to models. The measured val- 
ues are higher than those from the models considered, 
except for PYTHIA tune ATLAS-CSC for the 0.9 TeV 
and 2.36 TeV data, and PHOJET for the 0.9 TeV data, 
which are consistent with the data. At 7 TeV, the data 
are significantly higher than the values from the models 
considered, with the exception of PYTHIA tune ATLAS- 
CSC, for which the data are only two standard devia- 
tions higher. We have also studied the relative increase 
of pseudorapidity densities of charged particles (Table [TJ 
between the measurement at 0.9 TeV and the measure- 
ments at 2.36 TeV and 7 TeV. We observe an increase 



of 57.6 % ± OA%{stat.)tli %{syst.) between the 0.9 TeV 
and 7 TeV data, compared with an increase of 47.6 % ob- 
tained from the closest model, PYTHIA tune ATLAS- 
CSC (Fig. [T|). The 7 TeV data confirm the observation 
made in [HITS] that the measured multiplicity density in- 
creases with increasing energy significantly faster than in 
any of the models considered. 

In Fig. [21 we compare the centre-of-mass energy de- 
pendence of the pseudorapidity density of charged parti- 
cles for the INEL>0|^|<i class to the evolution for other 
event classes (inelastic and non-single-diffractive events), 
which have been measured at lower energies. Note that 
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Multiplicity A?^^ Multiplicity 



Fig. 3. Measured multiplicity distributions in |r;| < 1 for the INEL>0|^[<i event class. The error bars for data points represent 
statistical uncertainties, the shaded areas represent systematic uncertainties. Left: The data at the three energies are shown 
with the NBD fits (lines). Note that for the 2.36 TeV and 7 TeV data the distributions have been scaled for clarity by the 
factors indicated. Right: The data at 7 TeV are compared to models: PHOJET (solid line), PYTHIA tunes D6T (dashed line), 
ATLAS-CSC (dotted line) and Perugia-0 (dash-dotted line). In the lower part, the ratios between the measured values and 
model calculations are shown with the same convention. The shaded area represents the combined statistical and systematic 
uncertainties. 



INEL>0|^l<i values are higher than inelastic and non- 
single-diffractive values, as expected, because events with 
no charged particles in \ri\ < 1 are removed. 

The increase in multiplicity from 0.9 TeV to 2.36 TeV 
and 7 TeV was studied by measuring the multiplicity dis- 
tributions for the event class, INEL>0|^|<i (Fig. [3] left). 
Small wavy fluctuations are seen at multiplicities above 
25. While visually they may appear to be significant, one 
should note that the errors in the deconvoluted distribu- 
tion are correlated over a range comparable to the multi- 
plicity resolution and the uncertainty bands should be seen 
as one-standard-deviation envelopes of the deconvoluted 
distributions (see also [3]). The unfolded distributions at 
0.9 TeV and 2.36 TeV are described well by the Nega- 
tive Binomial Distribution (NBD). At 7 TeV, the NBD 
fit slightly underestimates the data at low multiplicities 
(A^ch < 5) and slightly overestimates the data at high 
multiplicities (TVch > 55). 

A comparison of the 7 TeV data with models (Fig. [3] 
right) shows that only the PYTHIA tune ATLAS-CSC is 
close to the data at high multiplicities (iVch > 25). How- 
ever, it does not reproduce the data in the intermediate 
multiplicity region (8 < N^h < 25). At low multiplicities, 
{Nch < 5), there is a large spread of values between dif- 
ferent models: PHOJET is the lowest and PYTHIA tune 
Perugia-0 the highest. 



Conclusion 

We have presented measurements of the pseudorapidity 
density and multiplicity distributions of primary charged 
particles produced in proton-proton collisions at the LHC, 
at a centre-of-mass energy ^/s = 7 TeV. The measured 
value of the pseudorapidity density at this energy is sig- 
nificantly higher than that obtained from current models, 
except for PYTHIA tune ATLAS-CSC. The increase of 
the pseudorapidity density with increasing centre-of-mass 
energies is significantly higher than that obtained with 
any of the used models and tunes. 

The shape of our measured multiplicity distribution is 
not reproduced by any of the event generators considered. 
The discrepancy does not appear to be concentrated in 
a single region of the distribution, and varies with the 
model. 
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